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ABSTRACT 



The explosion of sub-Chandrasekhar mass white dwarfs via the double detonation scenario is a potential explanation for type la 
supernovae. In this scenario, a surface detonation in a helium layer initiates a detonation in the underlying carbon/oxygen core leading 
to an explosion. For a given core mass, a lower bound has been determined on the mass of the helium shell required for dynamical 
burning during a helium flash, which is a necessary prerequisite for detonation. For a range of core and corresponding minimum 
helium shell masses, we investigate whether an assumed surface helium detonation is capable of triggering a subsequent detonation 
in the core even for this limiting case. We carried out hydrodynamic simulations on a co-expanding Eulerian grid in two dimensions 
assuming rotational symmetry. The detonations are propagated using the level-set approach and a simplified scheme for nuclear 
reactions that has been calibrated with a large nuclear network. The same network is used to determine detailed nucleosynthetic 
abundances in a post-processing step. Based on approximate detonation initiation criteria in the literature, we find that secondary core 
detonations are triggered for all of the simulated models, ranging in core mass from 0.810 up to 1.385 M Q with corresponding shell 
masses from 0. 126 down to 0.0035 M Q . This implies that, as soon as a detonation triggers in a helium shell covering a carbon/oxygen 
white dwarf, a subsequent core detonation is virtually inevitable. 

Key words, supernovae: general - nuclear reactions, nucleosynthesis, abundances - hydrodynamics - methods: numerical 



1. Introduction 

The basic physical mechanism for type la supernova (SN la) 
explosions has become widely accept ed ever since it was first 
propo sed almost 50 years ago (see e.g. lHillebrandt & Niemeveil 
2000, for a review ): a thermonuclear exp losion in electron- 
degenerate matter dHovle & Fowled [i960) produces radioac- 



. tive Ni that, b y its decay, release s the energy that p owers 
• ' the light curve (ITruran et all Il967t IColgate & McKeel 1 19691 
£^ | iKuchner et al.l fl994l) . Despite this long history, the questions 
pZ. • of the progenitor system and the explosion scenario have not 
^— ' been completely answered. In the so-called double detona- 
. . , tion sub-Chandrasekhar model, a detonation in an accreted he- 
k>" ■ lium shell causes a secondary deton ation of a carbon/oxygen 
white dwarf (C/O W D) core (e.g. IWooslev & Weaver! Il994t 
ILivne & Arnettlll995l) . This happens at a total mass below the 
Chandrasekhar limit and can lead to a wide range of possible 
explosion strengths. Recent population synthesis studies suggest 
that such events are in principle frequent enoug h to account for 
a significant part of the observed SN la rate dBelczvnski et al.l 
l2005llRuiter eTalEo ogfl 

This scenario hinges on two critical points - first the forma- 
tion of a detonation in the helium shell and second whether a 
successful detonation of the helium shell can detonate the core. 
Here we focus on the second question. A secondary detonation 
can be triggered in two different ways: either directly when the 
helium detonation shock hits the core/shell interface ("edge-lit"), 
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1 But note that there is still some discrepancy between these predic- 
tions a nd those made from observations, which tend to be significantly 
lower dDeloveetafll2b 05: Bildst en et alj2007h . 



or with some delay, after the shock has converged near the cen- 
ter. The edge-lit case is more restrictive since it requires a strong 
shock wave and might only work if the helium detonation starts 
at som e altitude above the base of the shell (cf. lLivne & Glasnerl 
1990). In this work we examine the delayed mechanism since 
it may still lead to a core detonation even for shocks too weak 
for the edge-lit case. In this spirit of determining minimal con- 
ditions for a core detonation, we look into models with helium 
shell masses that are substantially lower than previously consid- 
ered. Provided a detonation is triggered in these, the shock they 
dri ve into the core is expected to be particularly weak. 

iFink et alJ (120071 hereafter referred to as Paper I) found a 
secondary core detonation to be robustly triggered in multidi- 
mensional hydrodynamic simulations of generic 1-M models 
with shell masses of 0.1 to 0.2 M Q . Here we revisit the double- 
detonation sub-Chandrasekhar supernova scenario for a series of 
models with different core-shell mass combinations. We use the 
minimum heliu m shell mass e s req uired for dynamical runaway 
calculated by Bilds ten et al.l d2007l) . We find that even for this 
most conservative case, a helium detonation initiated at the base 
of the shell robustly triggers a secondary detonation in the C/O 
core. We therefore also investigate the hydrodynamic evolution, 
nucleosynthesis, and observational implications of such explo- 
sions. 

Section|2]contains the setup of our simulations. In Sect. [3] we 
describe the hydrodynamics and nuclear reaction network codes 
used to perform the simulations. The results are presented in de- 
tail in Sect. [4] Section [5] discusses the results and their obser- 
vational consequences in an astrophysical context. Our work is 
summarized in Sect. [6] 



1 



M. Fink et al.: Double-detonation sub-Chandrasekhar supernovae: can minimum helium shell masses detonate the core? 



Table 1. Parameters of our initial models. 



Model 


1 


2 


3 


4 


5 


6 


M tot 


0.936 


1.004 


1.080 


1.164 


1.293 


1.3885 




0.810 


0.920 


1.025 


1.125 


1.280 


1.3850 


M sh 


0.126 


0.084 


0.055 


0.039 


0.013 


0.0035 


Pel 


1.45 


2.4 


4.15 


7.9 


28 


140 


Pb,5 


3.7 


4.0 


4.5 


6.1 


6.4 


8.7 



Notes. p c 7 denotes the central density in units of 10 7 g cirr 3 , p b 5 is the 
density at the base of the helium shell in units of 10 5 g cm 4 . M tat , M core , 
and M sh are the masses of the WD, the C/O core, and the helium shell, 
respectively. All masses are given in solar masses. 



2. Explosion scenario 

2.1. Initial models 

The models we study are WDs with a composition of 12 C and 
ie O (equal parts by mass) with minimum accreted helium shell 
mass (M s h) that mig ht lead to a detonation, as calculated by 
Bilds ten et al.l (I2007I) . They assumed a certain M C0le -M s h com- 
bination and calculated the resulting structure and P^, T\, evo- 
lution (with Moore being the mass of the accretor WD, and Pb 
and 7b being the pressure and temperature at the base of the 
shell) in a series of hydrostatic equilibrium integrations assum- 
ing a constant-temperature core of r core = 3 x 10 7 K and a fully 
convective shell with an adiabatic temperature profile. Heating 
by the triple-a reaction in a pure helium shell was taken into ac- 
count. Thus, the maximum temperatures during the flash could 
be determined and also the shortest nuclear burning timescales. 
In series of calculations with varying M s h, the minimum flash 
ignition masses that allow for a dynamical burning event and 
potential detonation M m ; n were determined (the term dynamical 
refers to the condition that the nuclear burning timescale at the 
shell base, where it is hottest, becomes as short as the local dy- 
namical timescale). 

Bilds ten et alJ (T2007) also study AM Canum Venaticorum bi- 
naries consisting of a C/O WD accreting at high rates from a he- 
lium white dwarf donor. They show that for WD masses ^0.8 M 
flash masses reach or surpass M m \ n before the donor mass is de- 
pleted. This makes these systems interesting candidates for the 
sub-Chandrasekhar scenario as they might be frequent and also 
might explode at small helium shell masses (0.01-0.1 M ). 

The initial models for our hydro simulations were con- 
structed as follows: given p c (central density), r core , p\,, and 7b 
at the instant when burning becomes dynamical (provided by 
Bildsten and Shen, private communication), hydrostatic equi- 
librium was integrated u sing our own code's equation of state. 
It differs slightly from Bild sten et alj d2007l) . e.g. it does not 
include Coulomb corrections to the ion pressure. Therefore, 
our core/shell masses might slightly differ, but this should not 
change the main properties of the explosion dynamics. 

The most important parameters of our initial models are 
given in Table [TJ In our models we neglect any potential en- 
hancements of the chemical composition due to the metallicity 
of the progenitors. The core masses M core range from ~0.8 M 
for Model 1 to nearly the Chandrasekhar mass for Model 6. As 
already mentioned, the shell masses M S h are set to M m ; n . The 
density at the base of the shell varies only moderately between 
3.7 and 8.7 x 10 5 gcirT 3 . 



2.2. Detonation initiation 

The simulations are carried out in a similar manner as in the pre- 
vious study dPaper J) . In all models, the helium shell is ignited 
at a single point at its base. This one-point ignition scenario in- 
troduces an asymmetry which makes it harder to trigger a core 
detonation than in symmetric shell-ignited models. Moreover, it 
is expected to introduce viewing angle effects in the synthetic 
light curves and spectra. As the simulation is carried out in 2D 
rotational symmetry, the ignition spot is placed on the positive 
z-axis of the cylindrical coordinate system. 

We follow the detonations in the helium shell and in the WD 
core with two separate level set functions (see also Sect. I3.21 i. 
The helium detonation is started by hand by setting the first level 
set function to positive values in a volume of choice at the base 
of the shell. Since dynamical burning is a necessary but not suffi- 
cient criterion for the initiation of a detonation, the formation of 
the helium detonation is a fundamental assumption in our mod- 
els. 

After shell ignition, the WD core is scanned for "hot spots" 
arising from the shock exerted on the core by the helium det- 
onation. If a sufficiently large volume becomes hot and dense 
enough, a C/O detonation is initiated by setting the second level 
set function to positive values in this whole volumfl To de- 
cide when to initiate the detonation, critica l densities and tem- 
peratu r es given in the liter ature are used dNiemever & Wo oslev 
(ll997h : lR5pke et al.l (120071) : see Tables 1 and 2 of lPaper J) . 

These detonation initiation criteria, which are based on the 
assumption of spherically symmetric linear temperature gradi- 
ents extending into u nburned fuel, are limite d by some inher- 
ent uncertainties (e.g. ISeitenzahl et al.l [20091 point out that the 
functional form of temperature and density gradients of the hot 
spot significantly affects the critical detonation initiation condi- 
tions). Furthermore, the process depends on the composition of 
the fuel, the background temperature and local velocity fields^. 
Since the critical length scales pertinent to detonation initiation 
are generally much smaller than the grid resolution, addressing 
these details is currently impossible in our full-star simulations. 
However, the detonation initiations in this work are fairly ro- 
bust: the critical temperatures and densities we apply are signifi- 
cantly exceeded (see Table|2]i. Fully resolved calculations of the 
detonation initiation conditions are expected to only moderately 
change the critical conditions. Thus, our conclusions should hold 
despite these uncertainties. 

3. Numerical simulation 

Our numerical hydrodynamics scheme, brie fly summ arized be- 
low in Sect. 13.11 is similar to the one used in Paper I. However, 
there have been major changes to the treatment of the burning 
physics, which are described in Sect. 13. 21 Also, a post-processing 
scheme is now employed to determine more realistic nucleosyn- 
thetic yields (Sect.l331l. 

3.1. Hydrodynamics 

In our Eulerian hydrodynamics code the reactive Euler equa- 
tions are solved using a finite volume scheme based on 

2 To prevent grid geometry effects on the flame, every hot cell initi- 
ates a burning bubble with a radius of three cells around it. 

For such effects, volume burning behind the converging shock wave 
would have to be accounted for, but this is not yet implemented in our 
code. 
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the PROMETHEUS code by iFrvxell et all (Il989t) . which is 
an implementation of the "piecewise parabolic method" of 
IColella & Woodward! d!984l) . In order to track the expanding 
WD during the explosion, a co-e xpandin g uniform grid as in 
iRbpke & Hillebrandtl d2005l) . and iRopkd d2005h is usecfl In 
some of the simulations (Models 4-6) exponentially growing 
cell sizes (e.g. lReinecke et al.l2002bllRopke et all 2006) are used 
in the outer parts of the grid in order to keep a sufficiently high 
resolution in the C/O WD during the helium shell detonation and 
successive shell expansion. The grid resolution was 1024 x 2048 
cells in 2D rotational symmetry. An equation of state for WD 
matter and monopole gravity complete the system. The equa- 
tion of state includes contributions of an arbitrarily degener- 
ate and arbitrarily relativistic electron gas, nuclei that follow a 
Maxwell-Boltzmann distribution, photons, and electron/positron 
pairs. The dependency on changes in the number of electrons per 
baryon, Y e , is also taken into account. 

3.2. Nuclear burning 

The two detonations are followed by independent level sets in 
passive implementation ([Osher & Sethian| [l988; Reineck e et al.l 
1999; Golomb ek & Niemeverl 120051) . The burning velocity is 
calculated as a function of the local density: in the 4 He deto- 
nation the Chapman-Jouguet case is assumed and the burning 
speed relative to the ashes is the local sound speed. In C/O at 
high densities above ~2 x 10 7 self -sustained detonations are of 
pathological type (e.g. Sharpe 1999). For t he bur ning speed rela- 
tive to the fuel the values o flGamezo et al. ( 1999) are used above 
10 7 g cm 4 . At lower densities the Chapman-Jouguet case is as- 
sumed. 

A simplified scheme is used to mo del the energy release o f 
nuclear reactions, similar to lPaper J and lReinecke et all ((2002a): 
the code includes five species variables for 4 He, 12 C, ie O, inter- 
mediate mass elements (IMEs), and iron group elements (IGEs). 
Changes in composition and internal energy due to fast reac- 
tions occur instantaneously behind the flame discontinuity. To 
this end, the final abundances and g- va lues have been tabulated 
against fuel density. This new prescription constitutes an im- 
provement to the burning scheme of iPaperll Details are given 
in the Appendix. 

3.3. Post-processing 

Our post-processing scheme is similar to Travagl io et al. (2004): 
passively advected tracer particles are used in the hydro runs 
and detailed nucleosynthetic yields are calculated afterwards 
from their density and temperature trajectories. To this end, a 
large nuclear network with 384 species ranging from protons, 
neutrons, and a-particles to 98 Mo is employed. A description 
of the code used to solve the l arge nuclear network c an be 
found in Thieleman n et al.l d!990l) . lThielemann et al.ldl996l) . and 
Iwa moto et al.l dl999h . The reaction rate libraries adopted in this 
work are the same as in th ose references, however, newer tables 
for weak reaction rates dLanganke & Martinez-Pinedol [2000) 
were used as described in Travagli o et alJ ( 20041) . 

A number of n core = 80 x 160 = 12 800 tracer particles is 
placed equidistant in radial mass coordinate and cos 9 in the WD 



4 The expansion of a fixed mass shell is tracked. To resolve the shock 
convergence arising from the helium detonation, for Models 1-3 the 
grid is kept static until the end of helium burning and for models 4-6 it 
is kept static until the onset of C/O detonation. 



Table 2. Conditions at core detonation initiation for all models. 



Model 


tign [S] 




PS.Kn 


z len [km] 


|Zignl/^core 


A [km] 


1 


1.94 


4.62 


1.08 


-1910 


0.39 


9.86 


2 


1.78 


4.15 


1.58 


-1850 


0.40 


8.43 


3 


1.58 


4.42 


2.24 


-1750 


0.41 


6.81 


4 


1.30 


4.09 


1.89 


-1560 


0.41 


5.59 


5 


1.00 


4.14 


4.63 


-1370 


0.46 


3.86 


6 


0.64 


3.23 


8.92 


-1040 


0.50 


2.54 



Notes. fj gn , Tgj gn , psjgn, and z- Kn are the time, temperature in 10 9 K, 
density in 10 s g cirr 3 , and z coordinate at the ignition spot, respectively. 
i? core denotes the radius of the C/O- WD core and A the grid resolution, 
which is approximately 1/1000 of the WD radius. 



core such that each particle represents the same mas^l To spa- 
tially resolve the nucleosynthesis in the thin helium shell despite 
its lower mass, the same number of « s h = re core tracer particles 
was placed there. With this distribution, particles in the core and 
in the shell represent masses of -j 2 ^ and respectively. 

4. Simulation results 

In this section we present the main simulation results. All the 
models behave in a characteristically similar manner. Therefore 
our discussion will focus on Model 2 (M core = 0.92, M s h = 
0.084 M Q ) as a detailed example but we will also comment on 
the other models where appropriate. 

4.1. Explosive evolution 

Starting from the ignition point (see top-left panel of Fig. [TJ, 
the helium detonation surrounds the whole core until it con- 
verges again at the "south pole". While the detonation moves 
along the base of the shell, an oblique shock wave propagates 
into the c ore, converging off-cente r at a po int on the negative 
z-axis (cf. Livne & Glasner (1990); Paper ifl As can be seen 
in Fig. [T] this resembles the self-similar probl em of a spheri- 
cally or cylindrically conv e rging shock wav e (cf.lGuderlevl 19421 
Landau & Lifschitz 1991; Ponchaut et al. 2006), which results 
in strong shock compression if sufficiently small scales are re- 
solved. Due to self-similarity, the scales that need to be resolved 
to reach a certain compression just become smaller if the shock 
from the shell detonation is weake r. Howev er, the maximum pos- 
sible compression is limited (cf. Paper I). The question to be 
addressed in this study is whether the volume at which high 
enough temperatures and densities for dynamical carbon burn- 
ing are achieved is large enough for a detonation to form. In our 
simulations we always reached densities >10 8 gcirT 3 and tem- 
peratures >3.2 x 10 9 K on resolved scales >2.5 km. The critical 
radius for detonation formation at these conditions is, however, 
only 2 m. Thus, the critical conditions for detonation initiation 
are met for all our models, despite small shell masses. Table [2] 
lists t he cond itions at which the core detonations were ignited. 
As in iPaper fl. a conservative critical temperature of 4 x 10 9 K 

5 An offset is added to the coordinates such that each particle has a 
random initial position within its corresponding fluid element. 

6 Rotation effects that are neglected here could break the symmetry 
of the one-point ignition scenario if the ignition spot was located off the 
rotation axis. This could distort the minimum shocked volume and make 
a core detonation more difficult. The spherical shell ignition case and all 
ignition geometries that are symmetric with respect to the rotation axis 
should, however, not be influenced. 
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2 4 6 8 
r [1000 km] 



2 4 6 
r[1000 km] 



Fig. 1. Explosion evolution for Model 2. The 
density is color coded, and the solid cyan 
and magenta lines are the locations of the he- 
lium and C/O detonation flames, respectively. 
Dashed lines in cyan mark the border of the he- 
lium shell. 



was used. Thus, the given values are only a lower limit for the 
maximum possible compression at the given grid resolution A. If 
the detonation would have been suppressed, stronger compres- 
sion would have been achieved. Only Model 6 did not surpass 
4 x 10 9 K, despite being simulated at the highest spatial resolu- 
tion. The conditions reached in the shock compression, however, 
were still sufficient to safely assume a successful detonation trig- 
gering. Therefore, it was ignited at this lower temperature in a 
second run. This model verifies that it is harder to compress the 
core sufficiently if the helium shell mass is small. Conversely, 
the high initial density of the most massive model makes a det- 
onation easier. Based on our approximate initiation criteria we 
conclude that the limiting factor for a successful core detonation 
is only the successful formation of a detonation in the helium 
shell. 

Table [2] gives times and positions of the detonation initia- 
tions on the z-axis. That the carbon detonations occur earlier 
for smaller shell masses can be explained by the decrease of 
the core radii associated with the increasing core masses. At a 
smaller radius the helium detonation has a shorter way around 
the core while the helium detonation speed at the base of the 
shell is roughly constant for all models. Note that the approxi- 
mately self-similar nature of the problem results in the curious 



fact that the ignition spots of the core detonations lie at similar 
relative distances (0.4 - 0.5 ^ CO re) from the center. 



We will now discuss our fiducial case (Model 2) in detail. 
At densities $4 x 10 5 g cirT 3 in the helium shell, burning is rela- 
tively incomplete and nuclear statistical equilibrium (NSE) is not 
reached. The final composition is ~63 % of 4 He, ~ 1 0% of IMEs, 
and ~26% of IGEs (see Sect. 14.21 and Table [3] for more details). 

The C/O detonation starts at t ~ 1.8 s and at z 1900 km 

(see Fig. Q) and it produces 0.34 M of 56 Ni, 0.44 M G of IMEs 
and 0.11 M of le O. The structure of the ejecta at t = 10 s, 
where our simulations stop and the ejecta are close to homolo- 
gous expansion, is shown in Fig. [2] The shown distribution has 
been derived from the tracer particles after the post-processing 
step. Due to the two detonations, IGEs can be found both in the 
central region and in the shell. The shock of the C/O detonation 
partially penetrates into the helium-detonation ashes. This gives 
rise to the Richtmyer-Meshkov instability which generates vor- 
tices at the boundary between C/O- and helium-detonation ashes 
(see upper right part of Fig. [2]>. This effect causes some mixing 
between unburned C/O and IGEs from the helium detonation. 
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Table 3. Total nucleosynthetic yields of selected species or groups of species. 



Model 


1 


2 


3 


4 


5 


6 








C/O core detonation 








A^C/O.fuel 


8.1 x 10 


9.2 x 10 


1.03 


1.13 


1.28 


1.39 




Mice, 


1.8 x 10" 1 


3.6 x 10 _1 


5.7 x 10~' 


8.2 x 10-' 


1.11 


1.33 




MlMEs 


4.8 x 10" 1 


4.4 x 10 _1 


3.7 x 10~' 


2.6 x 10-' 


1.2 x 10~' 


3.1 x 10~ 


-2 




1.7 x 10~' 


3.4 x 10 _1 


5.5 x 10-' 


7.8 x 10-' 


1.05 


1.10 




M 5 2 Fe 


7.6 x 10~ 3 


9.9 x 10~ 3 


9.6 x 10~ 3 


7.9 x 10~ 3 


4.2 x 10~ 3 


1.7 x 10- 


-3 


M 48Cr 


3.9 x 10~ 4 


4.6 x 10~ 4 


4.5 x 10~ 4 


3.8 x 10- 4 


2.1 x 10~ 4 


7.1 x 10- 


-5 




1.4 x 1(T' 


1.1 x 10 _1 


8.0 x 10~ 2 


4.2 x 10~ 2 


3.1 x 10~ 2 


1.2 x 10- 


-2 


Ml2 c 


6.6 x 10~ 3 


4.4 x 10" 3 


2.7 x 10- 3 


8.8 x 10" 4 


5.9 x 10~ 3 


7.4 x 10 


-4 








Helium shell detonation 








He, fuel 


1.3 x 10~' 


8.4 x 10~ 2 


5.5 x 10~ 2 


3.9 x 10- 2 


1.3 x 10~ 2 


3.5 x 10- 


-3 


Ml GEs 


'i r\ - , 1 f\— 1 f i - * o frt \ 

2.9 x 10 ( 23%) 


2.2 x 10 ( 26%) 


1.7 x 10- 2 (30%) 


1.3 x 10~ 2 (33%) 


4.2 X 10 ( 32%) 


1.1 x 10" 


3 (31%) 




1 3 x 10~ 2 ( 10%^ 


8 ? x 10~ 3 C 10%^ 


5.3 x 10~ 3 ( 10%) 


5.7 x 10~ 3 ( 15%) 


1 9 x 10~ 3 C 14%"l 


7.3 x 10" 


- 4 (21%) 


M56 Ni 


8.4 x 10~ 4 ( 1%) 


1.1 x 10~ 3 ( 1%) 


1.7 x 10~ 3 ( 3%) 


4.4 x 10- 3 ( 11%) 


1.5 x 10~ 3 ( 11%) 


5.7 x 10" 


" 4 ( 16%) 


Afe Fe 


7.6 x 10~ 3 ( 6%) 


7.0 x 10~ 3 ( 8%) 


6.2 x 10~ 3 ( 11%) 


3.5 x 10~ 3 ( 9%) 


1.2 x 10" 3 ( 10%) 


2.0 x 10" 


" 4 ( 6%) 


M48 Cr 


1.1 x 10~ 2 ( 9%) 


7.8 x 10~ 3 ( 9%) 


4.4 x 10 -3 ( 8%) 


2.2 x 10" 3 ( 6%) 


6.8 x 10~ 4 ( 5%) 


1.5 x 10" 


" 4 ( 4%) 


M44 Ti 


7.9 x 10~ 3 ( 6%) 


5.4 x 10- 3 ( 6%) 


3.4 x 10" 3 ( 6%) 


1.8 x 10~ 3 ( 5%) 


4.9 x 10~ 4 ( 4%) 


6.2 x 10" 


" 5 ( 2%) 


M40 Ca 


4.7 x 10~ 3 ( 4%) 


3.2 x 10- 3 ( 4%) 


2.2 x 10~ 3 ( 4%) 


2.2 x 10~ 3 ( 6%) 


6.8 x 10~ 4 ( 5%) 


2.4 x 10" 


" 4 ( 7%) 


M* m 


8.4 x 10~ 2 ( 66%) 


5.3 x 10~ 2 ( 63%) 


3.3 x 10~ 2 ( 60%) 


2.0 x 10~ 2 ( 52%) 


6.9 x 10 -3 ( 54%) 


1.7 x 10 


3 (48%) 


M L 


2.0 x 10-' 


3.7 x 10- 1 


5.7 x 10' 1 


8.0 x 10- 1 


1.06 


1.10 




E km [10 51 erg 


] 0.90 


1.04 


1.20 


1.40 


1.59 


1.68 





Notes. Mc/o,fuci and Mnc.fuei are the total masses of initial fuel in the C/O core and the helium shell, respectively. For the helium detonation the 
values in brackets give the fraction of an isotope mass to the total shell mass M He _f ue i. M L is the total mass of all radioactive species that could 
power a light curve: 56 Ni, 52 Fe, and 48 Cr. All masses are given in units of M Q . is the asymptotic total kinetic energy. 



4.2. Nucleosynthesis 

The total nucleosynthetic yields of the explosions are given in 
Table [3] for each model. The upper part lists the results of the 
C/O core detonation, whereas the lower shows those for the he- 
lium shell detonation. Neutronization becomes important only 
for the highest mass model showing 16% of IGEs that are not 
56 Ni. This model is also peculiar in having almost no IMEs (only 
~2% of the total mass). The NSE freeze-out is a-rich for all 
models showing significant contributions of 4 He, 57 Ni, 58 Ni, and 
60 Zn in the final composition. Fig.[3]shows the distribution of the 
nucleosynthetic products in velocity space along the equatorial 
axis (this is representative of the mean for the whole explosion). 

A slice of Fig. [2] at the z = plane gives an overview of all 
nucleosynthesis taking place in Model 2 (Fig.|4|. Products of the 
C/O detonation are located below velocities of ~ 13 000 kms"'. 
The yields are: IGEs: 39%, IMEs: 48%, I6 0: 12% (see also 
Table [3]). For the low central densities present in most of the 
models, relatively large amounts of IMEs and 16 are found. 
However, there is almost no unburned 12 C. 

The helium detonation products are located at v J 
13 000 kms" 1 above th e C/O detonation layers. As discussed 
bv lBildsten et al.l d2007l) . the burning products at those low den- 
sities differ sign i ficant ly from previously published values (e.g. 
lKhokhlovlll984 1 1989b . For Model 2 the most abundant nuclei 
are unburned 4 He (63%), the IGEs 48 Cr (9%), 52 Fe (8%), and 
^Ti (6%), and the IMEs 40 Ca (4%), 36 Ar (4%), and 32 S (1%). 
The low-density helium burning regime is characterized by not 
reaching NSE and the fact that at the low maximum temper- 
atures a-captures are much faster than the triple-o' reactions. 
For the low initial densities in the outer shell, helium is mostly 
unburned. Deeper in the shell the higher densities increase the 
triple-a rate meaning that more 4 He is burned and higher max- 
imum temperatures are reached. Once a 12 C nucleus is formed, 
a-captures process it rapidly to higher mass numbers. Since the 



Coulomb barriers increase with mass number, this process stops 
at some maximum mass number depending on the local tempera- 
ture. This leads to an inwards increase of the mean mass number, 
peaking roughly at 48 Cr or 52 Fe. The final yields of Model 2 in 
velocity space (Figs. |4]and|5]) clearly demonstrate this trend. In 
the very innermost parts of the shell, where the initial densities 
are highest, the mass numbers are again lower. This is due to 
enrichment by carbon produced by triple-a reactions that take 
place before the onset of the detonation. This occurs because 
of the high temperatures at the base of the helium shells in our 
initial models. In a mixture of 12 C and 4 He, if the carbon mass 
fraction Xn c exceeds a value of ^, then there is not enough 4 He 
to form nuclei with mass number A or higher by a-captures in 
the course of a detonation passing through this matter. E.g., for 
Xn c > | one expects A max < 36. 

The initial triple-a burning in the hottest regions of the shell 
in the initial models provided by Bild sten et al.l d2007l) has not 
been taken into account in our hydrodynamic simulations, as 
they lack a nuclear network that could calculate this volume 
burning effect. It has, however, been considered in our post- 
processing step. This introduces some asymmetry in the results 
which reflects our choice of a pure helium shell composition and 
single spot ignition: while the detonation wave wraps around 
the WD, more and more 12 C is produced in the remaining shell 
material. For Model 2, the mass-fraction of initial I2 C which is 
reached at the equator by this volume burning is consistent with 
the m inimum mass fraction values given in IShen & Bildsten! 
d2009h . Above the equatorial plane there is less initial carbon, 
below there is more. Therefore our total yields should be roughly 
consistent with spherically symmetric or one-dimensional simu- 
la tions based on t heir in itial compositions. The values published 
in lBildsten et al.l (2007) from the one-dimensional detonation of 
a similar model, however, differ from ours as they assume a det- 
onation in pure 4 He. Thus it is not surprising that they reach 
56 Ni by a-captures and find it to be the most abundant burning 
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5 10 15 20 25 

v r [1 000 km/5] 



Fig. 2. Ejecta structure in velocity space for Model 2 at 10 s. 
Color coded is the mean mass number (averaged over mass: A = 
A/Xj, where A, and X/ are mass number and mass fraction of 
a nucleus i, respectively). 

product. The fact that they burn ~57% of the helium and we only 
about 40% (in Model 3, which is closest to their model) could be 
due to the different explosion dynamics of our one-point ignited 
and their whole-shell ignited detonation. 

4.3. Asymmetry effects 

In Fig. |2] asymmetries in both composition and ejecta velocities 
of Model 2 are visible. These are quantified in Fig. [6] which 
shows the abundances of the main groups of elements aver- 
aged for three different angular ranges: 0°^45° ("north"), 67.5°- 
112.5° ("equator"), and 135°-180° ("south"). The main differ- 
ences are: 

- Helium detonation products in the north extend over a wider 
velocity range than those in the south. As the ejecta are al- 




10 20 30 

v [1000 km/s] 



Fig. 3. Mass fractions of the main groups of nuclei in velocity 
space for all models at 10 s. Shown is the average of an an- 
gle range of ±n/32 around the equator. The velocity bin size is 
100 km s" 1 . 




Q F .... i .... i ......... . . . , : 

5 10 15 20 25 30 

v[1000km/s] 

Fig. 4. Mean atomic mass number in the equatorial plane (z = 0) 
for Model 2 at 10 s. Points give the values for individual tracer 
particles within 9 - n + n/32. 



ready homologous to very good accuracy at 10 s, this cannot 
change by stronger expansion of the southern ejecta. The ef- 
fect is also visible in the final distribution of tracer particles 
shown in Fig. [7] As every particle represents the same mass, 
it can also be seen that the southern regions are denser than 
the northern ones. This property might be important for spec- 
trum formation, as a wider range in velocity space increases 
the range of possible photon absorption frequencies. 
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Fig. 5. Distribution of the isotopes in the ejecta in velocity space 
for Model 2 at 10 s. Shown is the average of an angle range of 
+7r/32 around the equator. The velocity bin size was 100 km s _1 . 



0° - 45° z 




10 15 20 

v[1000 km/s] 



Fig. 6. Asymmetry of mass fractions in velocity space for 
Model 2 at 10 s. (Group) abundances are averaged for three dif- 
ferent polar angle ranges and over velocity bins of 250 km s _1 . 



- From north to south the amount of burned material de- 
creases. It is 48%, 37%, and 24%, for the northern, equato- 
rial, and southern directions of Fig. [6l respectively. Also the 
amount of synthesized IGEs is less (corresponding to lower 
maximum mass numbers in Fig.|2]i: it is 88%, 66%, and 51% 
of the burned matter, respectively. 

- At the south the C/O core has been compressed more 
strongly: the outer boundary of the core ejecta is found to be 
roughly at 14000 kms" 1 , 13000 kms" 1 , and 12000 kms" 1 
from the top panel to the bottom panel of Fig. [6] This is also 
consistent with the IGE abundance from the core detonation 
increasing from north to south (30%, 36%, 55%) at the ex- 
pense of IMEs (56%, 50%, 38%) and oxygen. 

- While the core ejecta extend to lower velocities in the south 
the IGEs reach significantly higher velocities there. 

- Especially in the northern hemisphere there seems to be sig- 
nificant mixing of unburned carbon and oxygen from the 



i i i i i i i 
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v, [1000 km/s] 



Fig. 7. Distribution of tracer particles in velocity space for 
Model 2 at 10 s. Shown are only particles which initially belong 
to the helium shell. 



core with IGEs from the shell detonation. This might be due 
to the Richtmyer-Meshkov instability (see Sect. 14.11 ). 

Qualitatively similar asymmetry effects appear for the other 
models of our study. 



*= 5. Discussion 



The primary goal of this work has been to investigate whether 
or not secondary core detonations can be tri ggered for scenario s 
with the minimum helium shell masses of Bildsten et al. (2007). 
We find that secondary core detonation conditions leading to a 
successful explosion of the WD are obtained in all of our six 
simulations. We now comment on the key observational impli- 
cations of our double-detonation models. 

Depending on the initial central density of the model, 
the C/O detonation produces nickel masses between 0.17 and 
1.1 M Q . In principle, this range is sufficiently broad to en- 
compass all major c lasses of type la supernovae (cf., e.g., 
IStritzinger et al.|[2006l for a sample of 56 Ni masses determined 
for 16 well-observed SNe la): the high mass end of sub- 
luminous SN 1991b^-like events (~0.07-0.17 M : Model 1), 
normal type la supernovae (~0.4—0.8 M Q : Models 2-4), and 
bright SN 1991 iff like explosions (-0.85-1.0 M Q : Models 4, 
5). In Fig.|8]we show angle-averaged ultraviolet-optical-infrared 
(UVOIR) light curves for all six of our models ( computed with 
the artis code; [Rromer & Sim 2009; S iml 120071 corresponding 
values are given in Table @). The light curves illustrate that not 
only does the peak magnitude vary significantly between the 
models, as expected from the differences in nickel mass, but 
that there is also significant diversity in both the rise times (see 
Table|4]i and the post-maximum light curve shape. This diversity 



7 iFilippenko etail \ 19921) : lLeibundgut etail d!993t) . 

8 iPhillips et all d!992h . 
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20 30 
Time since explosion (days) 

Fig. 8. UVOIR bolometric light curves for the models. Two light 
curves are shown for Model 1 : the first (solid line) shows the cal- 
culation for the complete model while the second (dashed line) 
shows the result obtained when the contributions of radioactive 
52 Fe and 48 Cr are neglected. 

Table 4. Rise times f pea k and peak absolute magnitudes M^"}^ for 
the UVOIR bolometric light curves of the models. 



Model 


1 


2 


3 


4 


5 


6 


fpeak [days] 


18.6 


18.6 


18.0 


15.5 


14.4 


13.9 


M bol 

peak 


-17.6 


-18.2 


-18.7 


-19.1 


-19.4 


-19.5 


Am*" 1 


0.64 


0.40 


0.50 


0.55 


0.55 


0.66 



Notes. 



is the change in bolometric magnitude between maximum 



light and 15 days thereafter. Due to Monte Carlo noise magnitudes are 
uncertain by ~0.1 mag. 



arises from the different distributions of the burning products in 
both the core and helium shell. In particular, there is a clear trend 
for faster rise times in the brighter models. This occurs since, in 
the brighter models, the 56 Ni-rich core material extends to higher 
velocities and the opacity of the outermost layers is less owing 
to the lower masses of the IGEs made in the helium shell detona- 
tion. The fainter models show a single, fairly broad UVOIR max- 
imum while the brighter models have an initial peak with a weak 
shoulder appearing around 30 days after maximum light. This 
shape is qualitatively similar to that obtained for UVOI R light 
curve s for standard SN la models such as W7 dNomoto et al.l 
11984 : see figure 7 of lKromer & Siml (120091) . Full details of our 
radiative transfer simulations, including complete sets of syn- 
thetic light curves and sp ectra will be presented in a companion 
study (iKromeretal .1120101) . 

For most of the models (1^1), the C/O detonation pro- 
duces a significant quantity of IMEs (~0.48-0.26 M ), as re- 
quired to account for the strong lines of e.g. silicon, sul- 
phur and calcium which characterize the maximum-light spec- 
tra of SNe la. The extreme Model 6, however, makes almost 
no IMEs and Model 5 yields only a rather small IME mass 
(~0.1 M Q ). This means that they are unlikely to be promis- 
ing candidates to account for real SNe la. However, they are 
still interesting as a demonstration that detonations of such 
small shell masses can still trigger a secondary detonation in a 
C/O- WD core. Note that possible initial com positions favoring 
more oxygen (and neon) for massive W Ds (cf.lDommgue z et alj 
l200U lGil -Pons & Garcia- Berrdl2001h are not considered here. 



Extrapolating from Seitenza hl et alJ (l2009h . we expect that core 
detonations would still be triggered for a composition of 30% 
carbon and 70% oxygen (in mass). For compositions with sig- 
nificantly lower carbon fraction like in oxygen/neon WDs deto- 
nation criteria are not available yet. Detonation conditions and 
nucleosynthesis in these stars should be investigated in future 
studies. The concerns regarding the WD initial composition are 
alleviated by the fact that our ignition spots are far above the 
center (see Sect. 14.11 and Table |2| where the carbon fraction is 
expected to be hig her than in the innermost region (see, e.g., 
lAlthaus et ai]|2005l) . 

The material produced in the helium shell detonation has 
important observable consequences. Although 56 Ni is not very 
abundant in our models, significant mass-fractions are predicted 
for the radioactive nuclei 52 Fe, 48 Cr, and 44 Ti. All of these 
could play a part in powering the supernova light curve (cf. also 
Bilds ten et al.l l2007. who consider single helium detonation su- 
pernova light curves). The yields of these nuclei are given in 
Tableland, for Model 1, are as large as 20% of the 56 Ni mass. 
The total mass of all 56 Ni, 52 Fe, and 48 Cr contributions is given 
as M L in Table [3] 

52 Fe and 48 Cr have relatively short decay times and release 
a similar amount of energy per decay as 56 Ni. Since they are 
located in the outer ejecta (see Fig. [3]), they can contribute to 
the early phase light curve of the models in which they are abun- 
dant. This is illustrated in Fig.[8]where, for Model 1, we compare 
the light curve computed including the energy released by 52 Fe 
and 48 Cr to that obtained if these decay chains are neglected. 
Generally, the contribution from 52 Fe and 48 Cr decays is fairly 
small and is most significant during the rising phase, as expected. 
The light curve at maximum is completely dominated by energy 
released from 56 Ni and 56 Co decays and it remains so through- 
out the decay phase. There is a modest enhancement around 
30 days (of about <0.2 mag) that is mostly due to energy pro- 
duced by 48 V, the daughter of 48 Cr for which the decay time is 
ri/2 = 16 days. The half-life of 44 Ti is too long for its decay 
to directly contribute to the early light curve. Nevertheless, the 
abundance of titanium is crucial for U - and B-band light curves 
and spectra since even small amounts contribute significantly to 
the opacity of the shell. 

The high velocity IGEs produced in the helium detonation 
impose an important constraint on the ability of our models to 
reproduce the early spectra of observed SNe la. SN 1991bg- 
like objects show titanium in their spectra, but our models are 
too bright to fit their light curves. Normal SNe la do not show 
clear signatures of IGEs at early epochs. Similarly, although 
there is evidence of IGEs affectin g the pre-maximum spectra of 
SN 1991T-like explosions (e.g. iMazzali et al J [19951) . these are 
inconsistent with our models since the important features there 
are associated with iron rather than the IGEs predicted to be 
abundant from our nucleosynthesis calculation^. However, to 
fully address the issue of whether the presence of heavy ele- 
ments created during minimum helium-shell detonations is in 
contradiction with observations requires detailed radiative trans- 
fer simulations and consideration of possible uncertainties in the 
helium-shell nucleosynthesis. In particular, we would like to em- 
phasize that the nucleosynthetic outcome of the helium detona- 
tion depend s on the enrichment of the shell with light a nuclei 
and 14 N (cf. lShen & Bildstenll2009l) . Thus, the resulting amounts 
of 52 Fe, 48 Cr, and 44 Ti could be significantly lower than found in 



9 Note that the decay-time of 52 Fe is so short, ~0.5 days, that it is 
expected to have almost completely decayed to 52 Cr before the light 
curve is bright. 
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this st udy. This will be discussed in a future study (Kromer et al. 
l2010h . 

6. Summary 

Motivated by the robustness of a second ary core detonation that 
was found in a previous study dPaper J) , the double-detonation 
sub-Chandrasekhar scenario for SNe la was re-investigated. This 
time we studied the case of minimum shell m asses that might 
detonate, as calculated bv lBildsten et al.l d2007l) - this is the most 
conservative case with regard to the question of whether sec- 
ondary core detonations will be triggered. 

In order to improve the accuracy of the explosion dynam- 
ics over that of Paper I, we switched to a more realistic detona- 
tion prescription in our numerical hydrodynamics simulations, 
including an energy release that was calibrated with a large nu- 
clear network. The same network was used in a post-processing 
step to calculate more detailed nucleosynthetic yields. 

Even for shell masses as low as 0.0035 M© and with more 
re alistic he lium detonation physics (not producing pure 55 Ni as 
in iPaper J) a successful core detonation was found to be very 
robust (for a discussion of uncertainties concerning the detona- 
tion criteria see Sects. | 2.2I and 14. U . Thus we predict that the 
"Ja" SNe of iBildsten et al.l (|2007) turn into brighter candidates 
for "proper" SNe la. Given the robustness of the core detonation 
and the fact that these events might be frequent enough to explain 
a significant part of the SN la rate, the double-detonation sub- 
Chandrasekhar scenario is a promising candidate for explaining 
some of these thermonuclear explosions. 
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Appendix A: Detonation tables 

In Paper t we assumed a certain transition density between burn- 
ing to NSE and incomplete burning in the C/O detonation. 
Additionally, burning was stopped below a roughly estimated 
fuel density. Therefore the total amount of IGEs and IMEs pro- 
duced was relatively uncertain. For this work we 'calibrated' our 
burning using a large nuclear network (see Sect. 13.31 ) in an iter- 
ation of explosion simulations and post-processing steps: we set 
up detonations spherically expanding from the center of a WD 
star and placed a large number of tracer particles along one axis. 
As starting point, the production of pure 56 Ni at all fuel densi- 
ties was assumed. The density and temperature profiles attained 
in this way were then used to determine the detailed nucleosyn- 
thetic yields in a post-processing step (see Sect. 13.31 ) with the 
large network. Thus the abundances of our five species (or the 
respective 2" va l ues ) could be tabulated against the initial den- 
sity of the unburned matter p u . This result was then used in an- 
other hydrodynamic simulation of the detonation and a second, 
more accurate abundance table was calculated from the post- 
processing. After several such iterations a self-consistent solu- 
tion was reached. The resulting table (see Fig. I A. 1 b is used in 
all the simulations presented in Sect. [4] and ensures a consis- 
tent energy release and molecular weight in the hydrodynamic 
simulations. Analogously a table for detonations in pure helium 
was determined (Fig. IA.21 . Although our calibration procedure 
does not account for potential changes in the detonation strength 




Pu [Q cm" 3 ] 

Fig. A.l. Mass fractions of the species in the C/O detonation ta- 
ble plotted against the density of the unburned fuel p u 




P u (9 cm" 3 ] 

Fig. A.2. Mass fractions of the species in the helium detonation 
table against p u . 12 C and ie O abundances are not shown, as their 
values are too close to zero. 

in realistic simulations, it still provides a reasonable approxi- 
mation. This was confirmed by post-processing our simulations 
presented in Sect. [4] which gave results consistent with the hy- 
drodynamic explosion simulations. 

In the C/O detonation (Fig. lA.ll i all major burning phases 
are visible: carbon burning, oxygen burning, and silicon burning. 
Above densities of about 2 x 10 7 g ctrT 3 NSE is reached. NSE is 
represented only by IGEs in this table. A better representation of 
NSE, a temperature and density dependent mixture of 56 Ni and 
4 He, is calculated in a different module of the code. Helium burn- 
ing (~Fig. IA.2t produces mainly IGEs and only a small amount of 
IMEs at the low densities in our shells (<10 6 gem 3 ). A high 
fraction of the initial helium is left unburned. 
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